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I. INTRODUCTION 



The outstanding problem in the theory of strong interactions is to understand how 
hadrons are built from quarks and gluons. The essential information on the internal structure 
of hadrons can be obtained from hard scattering processes. In such processes the duration 
of the interaction between a projectile and a fast moving hadron is small and the latter can 
be resolved into separate quarks and gluons (partons). Thus in hard scattering processes 
a hadron can be considered as a collection of non-interacting partons. This description of 
hadrons emerged from the study of deep inelastic scattering (DIS) and evolved into what is 
nowadays known as the parton model of hadrons and found its fundamental justification in 
QCD. 

The key ingredient of the parton model are the parton densities which encode information 
on the distribution of partons with respect to longitudinal momentum. The recent trend in 
this field is to study processes which can provide additional information (e.g. distribution 
of partons in the transverse plane, their angular momentum) about the internal structure of 
hadrons. Such information is encoded in so-called Generalized Parton Distributions (GPDs) 
which were introduced in Refs. in the analysis of Deeply Virtual Compton Scattering 

(DVCS). The GPD formalism is rather universal and applicable for the analysis of many 
different processes, ranging from completely inclusive to exclusive ones (for comprehensive 
reviews see Refs. p, la |7|). 

Since the GPDs are essentially nonperturbative quantities, no methods exist so far, which 
allow to obtain them directly from QCD. At present the only realistic way to get GPDs is to 
try to extract them from relevant experimental data Isl lol, llOI I (^y hich test only convolution 
of GPDs and specific kernels) and lattice data ll|, ll3|, [ij] (which only can provide 
moments of GPDs). Probably only global fits will be selective enough to really determine 
more than the most dominant GPDs. Such fits naturally rely heavily on Q^-evolution. This 
problem can be already treated within perturbative QCD. The equation governing the 
dependence of GPDs - the evolution equation - is well known. Since the GPDs are related 
to matrix elements of certain operators, the evolution equation for GPDs follows from the 
renormalization group (RG) equation for the latter. In special cases the GPD evolution 
equation reduces to the famous DGLAP [H H Q and ERBL [H 

evolution equations. 

However, the general structure of solutions of the evolution equation for GPDs was poorly 
understood so far. 



Several approaches [21|, |2j, |23|, were proposed in the past to get a solution for the 



general evolution equation for GPDs. All of them try to exploit the fact that Gegenbauer 
moments of GPDs have a simple scale dependence. However, the Gegenbauer polynomials 
do not form a complete set in the region where the GPDs are defined. This leads to the 
known problems when one tries to restore the GPDs by their Gegenbauer moments. The 
numerical algorithms based on these methods are rather cumbersome and ineffective (see 
for discussion Ref. j^). 

Recently, we presented a new method for solving the GPD evolution equation It is 
based on the correct incorporation of the symmetry properties of the evolution equation and 
provides a clear physical picture for the evolution. In Ref. [1] we considered the evolution 
equation for the gluon GPDs related to the matrix element of the twist two gluon operator 
and for simplicity neglected the effects of the mixing with quark-antiquark operators. In 
the present paper we give a detailed description of this approach for the example of the 
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flavor nonsinglet quark GPD and present the solution of the evolution equation for the 
flavor singlet GPD which takes into account mixing between quark and gluon operators. We 
would like to mention that an approach similar to ours was recently proposed in Ref. 0|. 
It will be briefly discussed in Sect. IIVI 

The paper is organized as follows. In the Sect. |n]we recall the deflnitions of GPDs and 
the corresponding evolution kernels. In Sect. IIIII the symmetry properties of the evolution 
equations are discussed. In Sect. |TV| we consider in detail the evolution equation for the 
isovector quark GPD and construct its solution. In Sect. |V|we present the solution of the 
evolution equation for the singlet quark and gluon GPDs. Sect. I VII contains the concluding 
remarks. In the Appendix some useful formulae are collected. 



II. BACKGROUND 

Throughout this paper we shall use the notations of Ref. j^. The GPDs are usually 
defined in terms of matrix elements of certain non-local light-cone operators. We shall 
consider the isovector operator 

0\zi,Z2) = q{zin)^+ry{z2n) (2.1) 

and quark and gluon isosinglet operators 

Q{zi,Z2) = q{zin)j^q{z2n) , (2.2a) 

g{z,,Z2) = G+^'^(^in)G'/'*(z2Ti). (2.2b) 

Here n is the light like vector, = 0, 7"*" = n ■ 7 and t"- {a = 1,2,3 ) are the Pauli matrices. 
In Eq. fj2.2ap a summation over flavor indices is implied. The real coordinates zi, Z2 specify 
the position of the operator on the light cone. As usual, we shall imply but do not display 
explicitly, the Wilson lines between the flelds at points zi and Z2- Taking matrix elements 
of these operators between hadron states one obtains GPDs related to these operators. For 
example, the pion isovector quark GPD is defined as follows 

ze^'^H{x,^,t) = I ^e^-'^''U^\p')m-z,z)\n%p)) , (2.3) 

where the kinematical variables are |5| 

P=^^ ^=^TT7^' i=^P-P) ■ (2-4) 

Henceforth we shall suppress the t— dependence of GPDs since it is irrelevant for the 
evolution. It is convenient to fix the normalization of the vector n by the condition 
(Pn) = = 1. With such a normalization the matrix elements of the quark and gluon 
operators, Eqs. ()2.1|1 . ()2.2a|l and (j2.2b|l . respectively, as well as the coordinates Zi, Z2 become 
dimensionless. 

The isosinglet GPDs are defined in a similar way (see e.g. Ref. |J] for the definitions). 
We shall consider the nucleon GPDs 

:F'ix, = / ^ e^^^^ {p'\ Qi-z, z) \p) , (2.5a) 
-^'(^>0 = j -e'''^p'\gi-z,z)\p), (2.5b) 
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The functions J^'^'^s) contain two different Lorentz structures and determine four scalar 
GPDs, H'^, and E'^,E^. As the Lorentz structure of GPDs is irrelevant for the dis- 
cussion of their evolution we shall treat the functions J^'^'^s) scalar functions. The GPDs 
H,J^'^,J-'^ have support in the x— interval [—1,1]. The skewness parameter ^ (Eq. (j2.4j) ) 
is also restricted to the interval [—1,1] by kinematic, (in DVCS processes C, > 0). It is 
standard to distinguish two different kinematical regions, |.^| < the so-called DGLAP 
region, and |^| > |x| - the ERBL or central region. In these regions GPDs describe different 
physical processes. In the central region a GPD describes the emission of a quark antiquark 
(gluon) pair from the initial hadron, while in the DGLAP region it describes the emission 
and absorption of a quark or an antiquark. 



A. Evolution equations: Isovector case 

All functions H, T'^, depend on the normalization scale /z at which the operators (j2.H) . 
fl2.2a|) and ()2.2b|) are defined. This dependence is governed by the renormalization group 
(RG) equation for the operators in question. We shall use the "coordinate space" version of 
this equation j2^. The RG equation for the operator O'^ at one loop order takes form 

+ /^(^?)^) ^"(^1' ^2) = \n (91(^1, ^2) , (2.6) 

where as is the strong coupling constant and the operator 7i, which we shall refer to as 
Hamiltonian, is an integral operator. Eq. ()2.6|) represents a concise description of renormal- 
ization of the local composite operators. Indeed, expanding the l.h.s. and r.h.s of Eq. (j2.6|) 
in a Taylor series in z^^z^ one reproduces the RG equation for the local operators. The 
Hamiltonian 7f (which encodes information on the mixing matrix) can be found in Ref. j26l |. 
It is convenient to represent it in the following form 

n = CpM, (2.7) 

m = -e^ + 1/2, (2.8) 

where Cp = {N^ — l)/2Nc and the integral operators H^,El2 are defined as follows 



1 a 



m_^]{zi,Z2) = / da-{2ip{zi,Z2) -ip{z"2,Z2) -(p{zi,z^^)) , (2.9a) 



a 

my,] {z,, Z2) = f da r d(3<f{z'^2, 4) • (2.9b) 
Jo Jo 

We use here the standard notations z"^. = Zia + Zka, a = 1 — a. As will be discussed later, 
both Hamiltonians are invariant under collinear conformal transformations and these two 
structures are the only structures which can appear in the Hamiltonian at one loop level (see 
Refs. ji^,!!^). Also we shall see in which way these terms are responsible for the different 
evolution effects. 

Let us now take the matrix element of both sides of Eq. ()2.6|) . Then we get an equation for 
the function ip°'{zi, Z2) = {out\0°'{zi, Z2)\in), which is an ordinary function (not an operator) 
of two real variables 2:1, 2:2- Obviously it satisfies the same RG equation, ()2.(jj) . independently 
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of the choice of in and out states. Of course, the properties of the function ip"'{zi, Z2) depend 
strongly on the choice of the in and out states. Taking the matrix element between vacuum 
state and one hadron state gives the hadronic wave function, while taking the forward matrix 
element results in the parton densities. Usually one has to know all of these functions (wave 
functions, usual and generalized parton densities) not in coordinate space but in momentum 
space. In momentum space the RG equations look different for different physical situations. 
But it is useful to remember that all of them originate from one and the same RG equation 
for composite operators. In some sense the only difference between the DGLAP-, ERBL- and 
GPD-type evolution equation lies in the initial conditions for the evolution. The functions 
which one wants to evolve have different properties in all these cases. Therefore we prefer to 
work in the coordinate representation and only in the last step perform the transformation 
to momentum space. 

For the isovector current we define the function ^i^{zi^ Z2) by 

^6'^V(;^i,^2) = (vr^(p')|(^'(^i,^2)k"(p)) . (2.10) 

It is related to the isovector GPD as follows 

.,(....)^2e-'«-«)/..e"<'-.H(.,0. (2.11) 

The subscript ^ of the function Lp indicates that this function has fixed total momentum. 
Under translation it transforms as ^^{zi + a, 2:2 + a) = e~'^^^"'(p^{zi, Z2). We also introduce a 
function given by the integral 

^^{z) = j dxe'^^'Hix^i) (2.12) 

Provided that the GPD is a smooth function of x the function $5(2;) (and ^^^{zi^ Z2)) vanishes 
fast for z = zi — Z2 ^ ±cxo. The function V2j(^i, 2^2), as a function of zi + Z2, is a plane wave. 
We shall present the solution of the evolution equation for the function H{x,^), however, in 
order to deal with well defined expressions on intermediate steps it is convenient to regard 
the convolution of (p^{zi, Z2) with some smooth function z/(^), 

-^v{Cm{zi,Z2) . (2.13) 

One can always assume that this new function ip{zi, Z2) vanishes sufficiently fast as zi^ Z2 
±cx). Eq. ()2.13p is a Fourier transform with respect to the total momentum ^, therefore the 
function ip{zi,Z2) satisfies the same evolution equation ()2.6|) . We postpone its solution till 
Sect. HYland discuss first the RG equation for the singlet quark and gluon GPDs. 

B. Evolution equations: Singlet case 

The quark and gluon operators ()2.2aj) . ()2.2b|) . mix under renormalization. However, the 
C— odd operator 

Q'{z^, Z2) = Q{zi, Z2) + Q{z2, zi) (2.14) 
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can not mix with gluons and evolves autonomously according to Eq. ()2.(i|l . contrary to the 
C— even operator which we choose as 

Q+(^i, Z2) = (Q(zi, Z2) - Q{Z2, z^)) , (2.15) 

which mixes with the gluon operator (j2.2b|l . To write down the RG equation in a compact 
form we introduce the vector notation O = {O'^, O^} where = Q and O'^ = Q~^. Then 
one gets 

^'§J; + ^(9)§A 0\z,,z^) = -^[K'^0%z,,z^). (2.16) 



The Hamiltonians TV"^ can be found in Ref. j26[. Let us notice that we use here the con- 
ventional definition of the nonlocal gluon operator which differs by a sign from that used in 
Ref. 13] • It results in a change of sign for the off-diagonal kernels. 
The Hamiltonian Ti^^ has the following form 

= (H? - H^) + QiV, + ^ • (2.17) 

The integral operators Hf , are defined as follows 

a" 

^{^]{Zi,Z2) = / da—i2(^{zi,Z2)-'^{z'^2^Z2)-'{){zi,Z2i)), 
Jo a 

where the weight function is uj{t) = 4(1 + 2r). 

The Hamiltonian W coincides with H, (see Eqs. (I2IZ|),(EH1)) 

Wi = n. (2.19) 

The off-diagonal Hamiltonians read jl^ 

n^'i = CpA-^ (I + 2 Hi) , (2.20) 

= njAMf, (2.21) 

where E[| is defined in Eq. ()2.9bj) and the Hamiltonian has the form ()2.18a|l with the 
weight function uj{t) = 1 + 3r. The operator A is the operator of multiplication by {zi — z^) 

\Ai^\{z^,Z2) = {z^ - Z2)v{zi,Z2) . (2.22) 
Further, in full analogy with the isovector case we introduce the functions 

fi{z,,Z2) = {AO\z,,Z2)\p). (2.23) 

They are related to the nucleon GPDs ()2.5aj) . (j2.5bj) as follows 

f9^zi,Z2)= e-*«("^+"2) J dxe'^^'^^T^ix,^), (2.24a) 
/|(2,,^2) = e-'«("^+"^^ y"c/xe""i2J^«'+(x,0, (2.24b) 
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where 2:12 = z\ — Z2 and 

T^^^{x,i)=%{T\x,i)-T\-x,i)) . (2.25) 

For later convenience we introduce the notations 

J-i(x,0 = T\x,i)-T\-x,i), (2.26a) 

^2(x,0 = T\x,i). (2.26b) 

Convoluting /| with a smooth function i^*(0 one gets the function 

f(^i,^2) = 1 '^y\iWz,.z,) (2.27) 

which vanishes for 2:1, 2:2 ±C)0. Let us note also that the function f^{zi, Z2) {/'^{zi, Z2)) is 
symmetric (antisymmetric) under permutation of variables, Zi ^ Z2 and this symmetry is 
preserved by the evolution. 



C. Local operators 

Let us discuss now the problems which arise in solving the evolution equation for GPDs. 
For the sake of simplicity we consider the isovector operator 0^{zi,Z2)- As we said above 
the nonlocal operator 0"'{zi,Z2) should be understood as the generating function for the 
local composite operators, 

00 

0%z,,Z2)= ^i^TOl^, (2.28) 

k,m=0 

where O'j^^ is [nD)^q{0)'~f^T°'{nD)"^q{0)/k\m\ and is the covariant derivative. Inserting 
this expansion into Eq. ()2.6|) one gets an infinite set of equations describing the renormal- 
ization of the local operators. Solving these equations one finds a set of the multiplicatively 
renormalized local operators of twist two. They are 

O%,,i0) = {tdX'-'m 7+r«cf g(0), (2.29) 

\d+ + d+J 

where C^^ is a Gegenbauer polynomial. These operators enjoy an autonomous scale depen- 
dence to one-loop order 

(9 (9 \ fy 

where 77V is the anomalous dimension of the operator. One can rearrange the expan- 
sion fTMi in the following way Q 

00 „i 
O'^iz^, Z2) = J2 BNi'^izi - Z2)'' / du{uuY+^ 0%{z^^2) ■ (2.31) 

N=0 
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Here Bn = 2{2N + 3)/{N + 1)! and = 0%^q{z). Inserting the expansion ^HH} into 

(P?Tn|l one gets 



^^{zu Z2) =v^e-^«(^^+^^) ^"^(^iV + 3) c^(0 ^n{zi2) , (2.32) 

Ar=o 

where Zi2 = zi — Z2 and 

^7v(^i2) = (e^i2)"'/Vjv+3/2(e2:i2) , (2.33a) 

cn{0 = I'/^Cf (^1^ H{x,0- (2.33b) 

Here Ju{z) is the Bessel function. The coefficient c^iO is related to the matrix element of 
the local operator 0% as follows 

^e^''ecN{0 = 2'''~\7r\p')\OU0)\7r^ip)) . (2.34) 

Taking into account the tensor structure of C^(0) one concludes that the coefficient ^^c^iO 
is a polynomial in ^ of degree (3]. This is the so-called polynomiality condition. Since 
the operator (9^(0), Eq. ()2.29p . is multiplicatively renormalized its matrix element depends 
on the scale in a simple manner 

c%'{0 = 4,'iOL-'''^'\ (2.35) 

where L = as{fi'i)/c(s{fJ'2), and 6o = ~ ^^f- Thus multiplying the coefficients cn{0 
Eq. ()2.32|) by corresponding exponents one gets the solution of the evolution equation ()2.6p . 



21 



The representation (j2.32p was obtained ffist in Ref. 

The series in ()2.32|) converges for any zu but it converges non-uniformly with respect 
to Zi2 what causes some problems. Namely, in the Fourier transform of (|2.32j) it is not 
possible to interchange the integration and summation. Indeed, let us note that the Fourier 
transform of the function "^nIz) 

r dze—^r^iz) = ^^C, o^ ^(^' -^")(^-^) Cf ( t) (2-36) 



(iv + i)(iv + 2) ^ ^ ^ V 

has support in the central region |x| < |^| only, while the Fourier transform of (f^{zi,Z2) 
is nonzero in the whole a;— region, [—1, 1]. In other words taking the Fourier transform of 
(j2.32p one obtains a formal representation for the GPD H{x,^) 

Hi^, = ^ f; ^iv c^(0 e{x' - e) (l - Cf (^l) , (2.37) 



N=0 

where 



1 2iV + 3 



2{N + 1){N + 2) ' 

but this series diverges if H{x,^) is nonzero in the DGLAP region |a;| > ^. Indeed, taking 
into account that for large and t > 1, C^'^{t) ~ (t + y/t"^ — 1)^ one concludes that the 
coefficients C7v(0 grow for large as and one gets a divergence. 
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FIG. 1: The coefficients log |f^| plotted as function of for n = 50. 



One can try to give sense to the series ()2.H7j] by interpreting it as a distribution, i.e. 
one has at first to carry out a convolution with some test function and then to take the 
sum provided that the latter converges. However, it is clear that the distribution defined in 
a such way, H{x,C,), is different from H{x,C,), since the former has support in the central 
region only. The problem with the Fourier transform of Eq. ()2.32j) is related to the fact that 
the Gegenbauer polynomials C^'^{x/^) are not orthogonal in the interval [—1,1] provided 
1^1 < 1. Nevertheless they form a complete set and the function H{x, ^) can be restored by its 
Gegenbauer moments CAr(^), (see Eq. ()2.33b|) ). The corresponding algorithm was suggested 
in Refs. [24, .SQJ and relies on the reexpansion of some system of orthogonal polynomials 
on the interval [—1,1], say Legendre polynomials, in terms of Gegenbauer polynomials, 

n 

Pn{x) = J2'-niO^''CT{x/0. 

N=0 

Obviously, one can also express the Legendre moments in terms of the Gegenbauer moments 

n 

Pn{0 = Y.''niO^''cr,{0- (2.38) 

Then the series ^^o(^"'" l/2)Pn(0-fn(a^) converges to the function H{x,^). Since the scale 
dependence of the coefficients ciy{^) is known (Eq. (j2.35|) ) one can calculate the coefficients 
PniOy a consequence the GPD H{x,C,), at any scale. However, an attempt to apply 

this scheme in practice reveals the following problem - the coefficients (£) grow fast with 
n. The explicit expression for the coefficients {^) can be found in |2ll. |24|. For the sake 
of clarity we consider the limiting case ^ = 0. In this limit one finds that C,^cn{C,) 2^cjvj 
where 

2r(A^ + 3/2) 



cn 



j dxx^H{x, 0) . 



The remaining dependence of the coefficients cn on N is then weak. The coefficients is 
nonzero only if and n have the same parity, { n — N = 2k ). Denoting = 2^ one 
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gets 

(-l)^(n + iV)! 



r 



2"A;!(n- A;)!r(Ar + 3/2) ' 
The logarithm of the coefficients |f^| as function of for n = 50 is plotted in Fig. Q One 



sees that in the sum 



N=0 

each term is large while the sum itself should be a number of order 0(1). Thus to calculate 
the coefficients pn with a reasonable accuracy one has to calculate the coefficients cat with 
"infinite" accuracy. It is clear that a numerical algorithm based on such a reexpansion 
sooner or later will have an accuracy problem. Thus, the methods for solving GPD evolution 
equations suggested in Refs. 21, 24, 1^^, although being mathematically correct, result in 



inefficient numerical algorithms. In what follows we suggest another method to solve the 
GPD evolution equation which relies heavily on the symmetry properties of the latter. 



III. SYMMETRY PROPERTIES 

As it is well known the classical QCD Lagrangian is invariant under conformal transfor- 
mations. This symmetry, however, does not survive in the full quantum theory due to the 
renormalization effects. Nevertheless, at one-loop level the counterterms inherit the sym- 
metry of the classical Lagrangian (for a review see e.g. Ref. j^). This means that if two 
local (gauge invariant) operators are related to each other by a symmetry transformation, 

= 6^0^, their one loop counterterms are related by the same symmetry transformation, 
AO"^ = 6^AO^. If one uses the formulation in terms of nonlocal operators this statement is 
translated into a statement about the invariance of the evolution kernel (Hamiltonian) with 
respect to the symmetry (conformal) transformations. 

We consider operators of a special type — the symmetric traceless operators of twist two 
(Eqs. (|2.ip . (|2.2j) ). This set of operators is closed under renormalization. The restriction 
of the full conformal group to this class of operators gives the so-called collinear confor- 
mal group which is the ordinary SL{2, R) group. The generators of the SL{2, R) group, 
S~^,S~,S^, satisfy the commutation relations 

[S+, S-] = 2S°, [5°, S^] = ±S^ (3.1) 

and can be chosen as 

S+ = z^d, + 2sz , S- = -9, , S° = zd, + s. (3.2) 

Here the parameter s, the conformal spin, specifies the representation of the SL{2, R) group. 
Under an infinitesimal SL{2, R) transformation a nonlocal operator transforms as follows 

5,0(zi, Z2) = e„ (5f + S!^) 0{z,, z^) , (3.3) 

where the spin operators should be taken in the representation corresponding to the confor- 
mal spin s = 1 for the quark operator (both singlet and nonsinglet), Eqs. ()2.1|) . ()2.2a|) . and 



10 



to the conformal spin s = 3/2 for the gluon operator, Eq. ()2.2bj) . For later convenience we 
introduce the following notations for the two-particle generators 



= {S'^ + S^)\s=i , (3.4a) 

= iS^ + S^)U=y,. (3.4b) 

Since, as was discussed above, the conformal symmetry survives at one loop level one gets 

[n,S^]=0 (3.5) 

for the Hamiltonian ()2.7|) . For the singlet case the Hamiltonian has the matrix form (see 
Eqs. (|2.16p ) therefore the symmetry relation takes form 

J2{n'',S^^] = 0, (3.6) 



fc=l,2 



where the matrix S is defined as 



s- 



(3.7) 



For the individual components Eq. (j3.6j) reads 

[W^, S^] = , [W, S^] = , (3.8a) 
W = W , H"^ = S;^ . (3.8b) 

To solve the evolution equations it is helpful to take this symmetry into account. 

Let us note that the transformations ()3.3|) are the infinitesimal form of the finite 
{SL{2,R)) transformations 

Oiz,, Z2) ^ (cz^ + d^^cz^ + d^'Oizl 4) , (3.9) 

where ^ = i^^^k + b)/{czk + d), ad — be = 1, and the spin s = 1(3/2) for the quark (gluon) 
operators. To proceed we recall some facts about the representations of the SL{2, R) group. 



A. SL{2,R) group 

The SL{2, R) group is the group of the real unimodular 2x2 matrices, 

9=(^l^^y ab-cd = l. (3.10) 

A detailed description of the representations of SL{2,R) group can be found in [s^, Issj ]. 
Since the SL{2,R) group is a noncompact group its unitary representations are infinite- 
dimensional. They can be organized into three series, namely the discrete series, the principal 
and supplementary continuous series ji^]. The last one will not appear in our analysis. 

The representation of the unitary principal continuous series is determined by two pa- 
rameters - the conformal spin s, which takes values s = 1/2 + ip, where p is a real number; 
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and by a discrete parameter, the so-called signature e = 0, 1/2. We need to consider the 
representations with signature e = only. In this case, the representation, which is denoted 
by TP, can be realized by unitary operators 

acting on the Hilbert space of functions of a real variable equipped with the scalar product 



{Hi^2)= / dzM^)Mz)- (3.12) 

J —oo 

The generators in this representation have the form ()3.2|) with conformal spin s = 1/2 + ip. 

The representations of the discrete series are characterized by the integer or half-integer 
conformal spin s, 5 = 1/2,1,3/2,... and by the transformation law 

T{g"')ij{z) = ^^—^{z') . (3.13) 

They fall into two classes usually denoted by Df. The difference between them lies in 
the space of functions they are defined on. The representations Df{D~) are defined on 
the functions of the real variable which can be continued analytically to the upper(lower) 
complex half-plane. The scalar product in these cases reads 

2s - 1 



TC 



{Hi^2)± = / d'z{2lmzy'-'M^) M^) ■ (3.14) 



'±lmz>0 

It is clear that functions from the representation space of Df can be represented in the form 
of a Fourier integral 

POO 

xP^{x)= dpe^'P^ijj^ip). (3.15) 



In the momentum representation the scalar product reads 



(^ilV^2)± = r(2s) / dpp'-''Mp)Mp)- (3.16) 

Jo 

One sees that in order to have a finite norm ip{p) has to vanish at p = 0. 

The tensor product of two unitary representations can always be decomposed into irre- 
ducible components. It will become clear later that we are interested in the tensor product 
decomposition of the representations of the discrete series, Df. This decomposition depends 
strongly on whether one multiplies the representations of the same type or not. Namely, the 
decomposition of the tensor product (or (g) D~) contains only representations 
of the same type 

oo 

Df^Df = J2(BDf^,^. (3.17) 

n=0 

In contrast, the tensor product Df ® Dj (or Dj ® Df ) can be decomposed into direct 
integral of the representations of the unitary principal continuous series 

/"OO 

Df®Df = dp® TP. (3.18) 
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Now let us try to understand the group theoretical properties of GPDs in more details. 
To be concrete we consider the quark GPD. In coordinate representation the GPD (f{zi, Z2) 
inherits the symmetry properties of the operator 0{zi,Z2), Eq. (|3.9|) . Since the conformal 
spin s = 1, this transformation corresponds to the tensor product of the discrete series 
representations of the SL{2, R) group. Let us represent the quark and antiquark fields by 

q{zn) = q'^ (zn) + (zn) , 

q{zn) = q^ (zn) + q^ (zn) , (3.19) 

where the components q~^ and g"*" contain the creation operators and thus, only positive 
Fourier harmonics, e**^^"-*^, p > 0, while the components q~ and q~ contain the annihilation 
operators (negative Fourier harmonics, e~*^^"^^, p > 0). Then, it is clear that e.g. the matrix 
element {p'\q~ (zin)^^ q^ {z2n)\p) as a function of two real variables, zi, Z2, admits an analytic 
continuation to the lower half-plane for each argument. Therefore, the GPD (f{zi,Z2) can 
be represented as the sum of four functions 

^{Zl,Z2) = (f~^'^{Zi,Z2) + <f^~iZi,Z2) + ^p~~^{Zi,Z2) + ^p"'{Zi,Z2) , (3.20) 

where (f"^ ~ {p'\q°'{zin)'y~^q^{z2n)\p) and a{(3) = ±. Each of the functions ip'^^lzi, Z2) trans- 
forms according to the tensor product of the representations D'^^D^, with s = 1. We notice 
here that the functions ip^^{zi, Z2) being transformed to the momentum representation have 
support in the central region, while the functions ^9^^(^1,2:2) in the DGLAP one. Since the 
form of the tensor product decomposition depends on whether a = j3 or a ^ P 
one has to decompose the functions 99"^ over different set of functions in these two cases. 

Let us note that our conclusion is in agreement with the solutions of the evolution equation 
in the two limiting cases, = 1 and ^ = 0. In the first one, ^ = 1, when the DGLAP region 
shrinks to zero, to solve the evolution equation one has to use the expansion in Gegenbauer 
polynomials, in agreement with the decomposition ()3.17|) . In the case of forward scattering, 
^ = 0, one should use a Mellin transform to solve the DGLAP evolution equation. This, in 
turn, agrees with the decomposition (j3.18p . 

All this shows that the two kinematical regions, "DGLAP" (|x| > |^|) and "central" 
(I a; I < 1^1) are quite different, both from the physical and mathematical point of view and, 
therefore, should be treated in a different way. 



IV. QUARK GPD 

In this section we consider in detail the case of the quark GPD (j2.3p . We shall construct 
the solution of the evolution equation for H{x,C,)- As it was explained earlier it is convenient 
to start analysis from the evolution equation in the coordinate representation 

(^^,-^ + P{g)^^ ^{z,,Z2) = [nv]{z^,Z2) , (4.1) 

where the Hamiltonian 7i and function (f{zi, Z2) are defined in Sect. Ill Al 

We remind that the conventional strategy to solve equations of this type is to expand the 
function ip{zi,Z2) in eigenfunctions of the Hamiltonian H, ip{zi,Z2) ~ YliE^Ei^E- If such an 
expansion is found, the dependence of the expansion coefficients ce on the scale can be 
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easily restored and, provided that the corresponding sum can be effectively calculated, this 
determines (/? at any scale /i. We shall try to follow this scheme. 

We start the discussion with the following remark. In order to talk rigorously about 
expansion in eigenfunctions of the Hamiltonian Ti, one should specify the Hilbert space in 
which the eigenvalue problem for the operator has to be solved. In other words one has to 
define the scalar product on the space of functions ip{zi,Z2)- It is clear that this construction 
(scalar product, Hilbert space) is external with respect to the problem. Indeed, Eq. ()4.1|) 
is an integro-differential equation which knows nothing about the Hilbert space. However, 
to solve this equation it is convenient to introduce the structure of a Hilbert space, i.e. a 
scalar product, and its choice is almost completely in our hands. 

The only condition which restricts the choice of the scalar product is the requirement that 
the function (p{zi, z^) belongs to the Hilbert space, i.e. it has to be normalizable with respect 
to the chosen scalar product. This is the place where physics imposes some constraints. The 
other requirements (such as the SL[2, R) invariance in the case under consideration) are 
useful but not mandatory. 

We choose the following SL{2, R) invariant scalar product 



{ipi\ip2) = J dZidZ2{zi - Z2) (pi{Zi,Z2)(p2{Zl,Z2) . (4.2) 

Obviously, it is invariant under SL{2, R) transformations 

M^i, Z2) (cz^ + d)-^'{cz2 + d)-^'ipkiz[, 4) , (4.3) 
with s = 1 and z' = {az + b)/ {cz + d). Going to the momentum representation 

ip{z^, Z2) = J ^e'(-«(^^+^2)+"(^i-^^»(^(a;, (4.4) 

one gets for the norm of the function (f 

\M' = l j didx\dMx,0?- (4.5) 

Then, taking into account Eqs. (j2.11|) . (|2.13|) one concludes that the physical GPD has 
a finite norm ()4.5|) if the integral j dx\dxH{x,^)\'^ is finite. Since the GPD H{x,C,), as a 
function of x, has support on the interval [—1, 1] this implies that it should vanish faster 
than (1 ±x)^/^ at the endpoints and be continuous inside this interval. These requirements 
are in agreement with the standard assumptions on the properties of quark GPDs (see e.g. 
Refs. HEJ). 



Let us remark here that the other possible choice for the invariant scalar product 



1^(^,01 



2 



<^ll ~ / d^dx '^l 'V , > (4-6) 



is ruled out. Indeed, in this case the function ip{x,C,) in order to be normalizable has to 
vanish at x = what is physically unacceptable. 

Next, it is useful to introduce a new function ip{zi, Z2) related to ^p{zi^ Z2) by the simple 
relation 

i^{zi, Z2) = [A(p] {zi, Z2) = {zi - Z2)ip{zi, Z2) . (4.7) 
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One can easily check that if the function Lp transforms according to ()4.3j) the function i\} 
obeys the same transformation law with s = 1/2. The map A : (f ^ ip is a one to one 
isometric map of the Hilbert space defined by the scalar product (|4.2|) to the standard, 
L^(i? X R) Hilbert space, 



(^i|^2) = J dzidz2 ^i{zi, Z2) iJ2{zi, Z2) . (4.8) 

The transformation of the Hilbert space L^(i? x R) 

^(2:1, Z2) (czi + dy\cz2 + d)-^^{z[, 4) , (4.9) 

determines the unitary representation of the SL{2, R) group which is equivalent to the tensor 
product of the representations 



1 1 



According to the Eqs. ()3.17|) and ()3.18p the decomposition of this representation into irre- 
ducible ones has form 

pi^ = V + 2 / dp® TP. (4.11) 

n=0 



A. Eigenvalue problem 

We remind that our purpose is to construct the expansion of the GPD (p{zi, Z2) in terms of 
eigenfunctions of the Hamiltonian ()2.7j) . Unfortunately, one can check that the Hamiltonian 
TC (Eq. ()2.7|) ) is non-hermitian with respect to the scalar product ()4.2p . It is possible to 
find its eigenfunctions, (it will be done later) but they do not form a complete set and, 
in general, it is not clear how to construct the expansion over this set. Nevertheless such 
an expansion can be constructed. To do so we use the following trick. As was explained 
before, the Hamiltonian commutes with the generators of the SL{2, R) transformations (see 
Eq. ()3.5p ). and as a consequence it commutes with the Casimir operator of the SL{2,R) 
group, [n, J2] = 0, 

J2 = (s, + 82)' = 5+5r2 + SUS'u - I) , (4.12) 

where ^'"2 = 5" + 6*2 and / is the identity operator. Taking into account the explicit form of 
the generators ()3.2p one finds that in the case si = S2 = s the Casimir operator takes form 

J2 = -4-2^ 5,52 4% (4.13) 

where 2:12 = Zi — Z2- For s = 1 this operator is a self-adjoint operator with respect to the 
scalar product ()4.2|) (and for s = 1/2 with respect to the scalar product ()4.8|) ). So one can 
expand any function over the set of eigenfunctions of the Casimir operator. Nevertheless, 
the operators li. and J^, despite the fact that they commute, have different eigenfunctions. 
However, one of the crucial points of our approach is that, as we shall show, the expan- 
sion in eigenfunctions of the Casimir operator can be transformed into an expansion in 
eigenfunctions of the Hamiltonian. 
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So as a first step we find the eigenfunctions of tlie Casimir operator fl4.1Hj) for tlie con- 
formal spin s = 1/2. This problem is equivalent to the decomposition of the representation 
of the SL{2, R) group into irreducible ones and its solution is well known [s^. The Casimir 
operator has eigenfunctions of both the discrete and continuous spectrum. The eigenvalues 
of the Casimir operators are usually written in the form j{j — 1). The parameter j is called 
conformal spin and takes integer values for the eigenfunctions of the discrete spectrum and 
is j = 1/2 + ip, with p being real and positive, for the eigenfunctions of the continuous 
spectrum. 

For our purpose it is convenient to choose the eigenfunctions in the form iIj{zi,Z2) ~ 
g-«C(2i+22)v]>(^2;i2). Inserting this function into the equation 

3^^{z,,Z2)=j{j-l)^{zuZ2) (4.14) 

one finds that the function 2:^2^ ^^(^12) satisfies the equation for Bessel functions, 
J±ij-i/2){\(,\zi2)- It is therefore convenient to define the function 

^(z) = e-i(^-i/2) z-'/' J,_i/2 mz) . (4.15) 

For integer j it is a single valued function of z in the whole complex plane. For non-integer 
j it has a cut from to —00. 

The eigenfunctions of the discrete spectrum, Pj{zi, Z2), j = n + 1, have the form 

if (^1, Z2) = e-*«(^i+^^)vl/|(^i2) . (4.16) 
Calculating the scalar product one finds 

(/,eU0 = 27r5(e-O^, (4.17) 

where we used the standard notation for the scalar product of and Py . The subspace 

spanned by the function Pj(zi,Z2), ^ > (^ < 0) corresponds to the subspace D~ {D~^) in 
the tensor product decomposition ()4.11|) . 

The eigenfunctions of the continuous spectrum, P^'^^zi, Z2), j = 1/2 + ip, (there are two 
eigenfunctions corresponding to each value of j) have the form 

P/'±(^i, Z2) = e"*«(^^+^^)vl/|'^(zi2) . (4.18) 
The functions \E'|'^(-2i2) are defined as follows 



¥^'^{z) = 
■' 2 cos Tcj 



¥^{z^)-¥,Az^)] , (4.19) 



where z± = ±z + iO. Again, the subspaces spanned by the eigenfunctions Pj'^{zi, Z2) 
correspond to two copies of the subspaces in the tensor product decomposition (|4.1H) . 
We notice also that these two functions differ only by permutation of their arguments, 

Pf^^{z,,Z2)=Pf'-{z2,Z,). 
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They are also invariant under the interchange j I — j, Pj'^{zi, Z2) = Pfr^(zi, Z2). Taking 
into account that the function is proportional to a Hankel function 

one can check that Pj'^{zi, Z2) is an analytic function of Zi in the upper and of Z2 in the lower 
complex half-plane. Next we conclude that from the properties of Bessel functions follows 
that the functions ^>^'^[z), as functions of j', are entire functions in the whole complex plane. 
For a fixed j they vanish like 1^1"^ for \z\ —>■ 00. 

Calculating the scalar product of the eigenfunctions of the continuous spectrum one gets 

= n6^a'm-mip-p')+Sip + p')]^coth7rp, (4.20) 

where a (a') = ± and j = 1/2 + ip. 

The eigenfunctions fl4.16p . ()4.18|) of the Casimir operator form a complete basis of the 
Hilbert space L^(i? x R). Thus any function from L^(i? x R) can be expanded as follows 

-i 5^^(j>5(j)if (^1,^2) - t / djuj'^{j)af{j)P^'^iz„Z2) . (4.21) 

-00 Jl-ioo J 

Here uj{j) = 2j — 1, u;^{j) = (j — 1/2) cotyrj and the expansion coefficients a^(j) and a^{j) 
are given by the scalar products 

adj) = ij^m , = {J,^,±m ■ (4.22) 

It follows from the properties of the eigenfunctions that the coefficients af{j) are entire 
functions of the conformal spin j in the whole complex plane satisfying the symmetry relation 
af{j)=af{l-j). 

Using the explicit expressions for the eigenfunctions ()4.16|) . ()4.18|) one can represent the 
expansion coefficients ()4.22|) in the form 



afU) = X j^^^ dxP,.i (^^^ i^{±x,0 , (4.23) 

where Pj{x) are the Legendre functions and 

X = ie-Z^v^V^. (4.24) 

ip{x,C,) is the Fourier transform of ipi^iy ^2) as in 1)2.111) . It can be shown that if the function 
ipix^^) behaves like (1 — x)° for x — 1, then the coefficients a^{j) vanish as j-^/a-o ^]-^g^ 
j goes to infinity along the imaginary axis. 

In the following we get rid of the Fourier integral over ^ and construct the expansion for 
the function '^/'^(-z) which is defined by 

^(^1, Z2) = I ^e-«^i+^^) ^^zi - Z2) . (4.25) 
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Taking into account the symmetry property of the expansion coefficients, a^{j) = a^(l — j), 
and using Eq. ()4.19p one derives 

^^(z) = ^^(j>5(j)^|(z)-^ r ^uj{j)af{j)¥^{z^). (4.26) 

Let us discuss the properties of the expansion ()4.2(j|l . First, we note that when j goes to 
infinity along the imaginary axis the functions '$^j{z±) behave as 

, p-if 0-1/2) 

VI/ ~ 2-'/'——-— (4) (1 + 0{l/j)) (4.27) 



r{j + i/2)\2J 

so that the integrals in (j4.26p converge. Moreover, shifting the integration contour to the 
right we further improve the convergency of the integrals. Second, let us remind that the 
function ip under consideration has the form ()4.7p . ip^{z) = zip^{z), and therefore vanishes 
at ^; = 0. For small z, ^'l(-z) ~ z^~^, thus all terms in the sum except the ffist one with 
j = 1 vanish as z — > 0. The integral over j follows the line Re j = 1/2. In order to study 
its small z behavior let us shift the contour of integration to the line Rej = 3/2. Since the 
integrand has a pole at j = 1 due to sinvrj one should evaluate the residue at this point. 
Adding the contributions from residues and the ffist term in the sum one gets 

h(l) + a+(l) + a^-(l)] ^=,{z). 

The term in the square brackets is zero, what is easy to see using Eq. (j4.23p and taking 
into account that ips.i^) = idx^^{x). Thus one gets the following expansion for the function 

= ij^{z)/z 

ni^) = |g ^(j>^(j)^|(z) - ^-^^0>«^0')^|(^±)} ' (4-28) 

where the integration follows the line Rej = k, 1 < k < 2. We shall show below that this 
is an expansion in eigenf unctions of the Hamiltonian. But before we want to notice that 
due to the asymptotics ()4.27j) the integration contour in ()4.28|1 can be closed in the right 
half-plane. Calculating the integrals in ()4.28|) by residues one finds 

00 

^^{z) = z''J2^(^)Mj)'^%^)^ (4-29) 

where 

= ^^{-ly f'^Cf_l(^]Mx). (4.30) 



The representation ()4.29|) coincides with ()2.32|) . which was obtained in Ref. 2l|. Next, the 
sum in fl4.28|) can be rewritten as an integral over j resulting in a representation similar to 
that derived in Ref. 
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FIG. 2: The integration contour Ca in Eq. (|OT|l . 

where _^ 

= 1^^^ rfxP,_i (^^^ [5x^f](±x) . (4.32) 

The integration contour Ca is shown in Fig|21 Deriving Eq. ()4.31|) we represented the function 
as sum of two functions ^~^{x) and ip~^{x), = ip'^{x)+ip'^{x), such that f'^{x) = 

at a; <— 1^1, and = at x > |^|. We also notice that the contour can not be deformed 

to be parallel to the imaginary axis, since in this case the integral over j starts to diverge. 
All three expansions ()4.28j) . ()4.29j) . ()4.31|1 are equivalent, the difference appears when one 
wants to make a Fourier transform. One can interchange the Fourier transformation with 
integration (summation) over j only in the representation (j4.28j) . 

Let us now show that the expansion ()4.28|) runs over the eigenf unctions of the Hamilto- 
nian. First of all, we note that the Hamiltonian (j2.7p is not a hermitian operator with respect 
to the scalar product ()4.2|) . Then the commutativity, [H, J^] = 0, does not imply that the 
operators share the same eigenf unctions. Indeed, one can easily check that the eigenfunc- 
tions of the continuous spectrum (j4.18j) of the Casimir operator are not the eigenfunctions 
of the Hamiltonian ()2.7|) . 

Nevertheless, one can find functions which diagonalize the Hamiltonian. (Of course, they 
are not mutually orthogonal and do not form a basis of the Hilbert space.) To solve the 
equation 

JVj(^i, Z2) = -did2z'^2^j{zi, Z2) = j{j - l)^jizi, Z2) 

we use the ansatz (pj{zi,Z2) = e"*^*-^^^^2^<l'|(2;i2), where j is an arbitrary complex number. 
This results in a second order differential equation for the function $|(-2). It has two inde- 
pendent solutions $|(-2) = z"^\&|(z) and ^l_j{z) = 2~^^f_^(^). Due to commutativity of 
the integral operator HI and the differential operator one concludes that 

[U^j]{zi, Z2) = A{j) ^j{zi, Z2) + 5(j) <^i_j(2;i, Z2) , (4.33) 

whenever the action of the Hamiltonian HI on the function (pj{zi, Z2) is well defined. Substi- 
tuting 'Pj{zi^ Z2) in the form e~*^''^^"'"^2''$|(zi2) into (|2.9p one finds that the integrals converge 
when Re j > 1. Next, let us notice that when the variables a, (3 run over the integration 
region, the argument of the function <l>| varies from to z. Thus to fix the coefficients A{j) 
and B{j) it is sufficient to study the asymptotics 2 of the r.h.s and l.h.s. of ()4.33p . 
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In this case one can substitute ^j{z) by its leading term ~ ^ and gets B{j) = and 
A{j) = E{j) with 

EU) = 2 im - ^(2)] - -T-^T + I . (4.34) 

Furthermore, we have to fix the ambiguity in the definition of the function $|(zi2) related to 
the existence of a branching point at Zi2 = 0. To this end we note that the Hamiltonian ()4.33|) 
transforms the functions with support in the region 2:12 > (212 < 0) into functions with the 
same support property. Thus, we conclude, that both functions, 6'(±2;i2)e~*^^^^"'"^^^$|(±^i2), 
(or their linear combinations such as, e~*^^^i"'"^2^$|(±2;i2 + iO)) are eigenfunctions of the 
Hamiltonian corresponding to the same eigenvalues, E{j). We notice also that since the 
integral / dzz'^\'^^j{±z)\'^ is finite for Re j > 1/2, the eigenvalue E{j) belongs to the discrete 
spectrum of the operator H. 

Thus we have shown that for any complex j such that Re j > 1 the function 

iffiz,, Z2) = e-'^^'^+''h^i¥^i±z,2 + tO) (4.35) 

is an eigenfunction of the discrete spectrum of the Hamiltonian EI with the eigenvalue given 
by dOl. 



B. Solution of the evolution equation 

Since we have shown that the expansion in the Eq. ()4.28j) contains eigenfunctions of the 
Hamiltonian one can easily write down the solution to the evolution equation ()2.6|) in the 
LO approximation 



[j=2 

[^^^^ ^^U)arU)L-^^'^^{zA , (4.36) 
^2j3_.^ smvrj « ^ | 



where a;(j) = 2j — 1, 



"^^1 and 7(j)=2^i?(j)- (4.37) 



a(/i') 60 
After Fourier transform Eq. ()4.36|) can be cast into the form 

(X) = E <0-) pS" {£) + E « /J ^'"(i)^-""?^ (f ) . (4.38) 

where the integration follows the line Rej = fc, 1 < k < 2. The expansion coefficients are 
given by the following expressions 

<(j) = V, 1'^^^'^ dxP,., d.^l{x) , (4.39a) 



c'rU) = dxP.^i (^^ ) [d.vl]i±^) , (4.39b) 



20 



where Pj-i{x) are the Legendre functions of the first kind and vj = {2j — l)/2. Notice, 



that the coefficients Ct'^(j) are antisymmetric under the interchange j — > 1 — j, Ct'^(j) 



— c^'^(l — j). We also want to stress here that the coefficients c^'^(j) are entire functions of 
j in the whole complex plane. 

The functions pf\x), qf\x) are expressed in terms of the Legendre functions of the first 

aud second kind We give here the expre^ions for tire furretions '(x), ,<™'(.) for 
general integer m since they appear in the solution of the evolution equation in the singlet 
sector. These functions (up to some normalization factors) are the Fourier transform of the 
functions z-^'^-iz) ipf'\x)) and z-^'^-iz + iO) The function pf^\x) is defined 

as follows 



pf\x) = {-ire{l - x^){l - x'r^'Pplix) = ry{l - x'){l - x'rc;'_V,L\{x) , (4.40) 

r(2m + 1) r(j - m) 



^3 

where 



' ^ ' r(m + l) r(j+m)' 



The function qj^\x) = for x < — 1, while for x > — 1 it is given by 

g(™)(a;) = (x^ - ir^^Q-^Hix), X > 1 (4.41) 



g(-)(x) = — ,{i-x'r/^Pr-^{-x), \x\<l. (4.42) 

■' 2 sm TTj 



The function qf\x) is continuous at the point x = 1, 



^f(-)U=-JJ^y (4-43) 
but its first derivative has a logarithmic singularity at this point 

-^gf^(l + x) ~ -^logx + ... . (4.44) 

ax ■' 2 

The formula ()4.H8j) represents the solution to the evolution equation for the pion isovector 
quark GPD, H{x,C,) = It can be used both for numerical and analytical study of the 

evolution. 

Let us discuss the solution ()4.38p in more details. First of all we notice that the integrals 
over i in (j4.38|) vanish whenever |x| > 1. Second, at the input scale /x' = /i (L = 1) the 
Eq. ()4.38|) can be represented as an expansion in eigenfunctions of the Casimir operator. To 
this end we shift the integration contour to the line Rej = 1/2. Since the function qf\x) 
for |x| < 1 has a pole at j = 1 one should calculate the residue at the point j = 1. Taking 
into account that 

1 



(j)l.=i = T-V^^(±|^|) (4.45) 



and (1 — x^)^/^Pn ^(""3;) = (1 + x) one finds that the contribution from residues is 

(4.46) 





X 
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Using the antisymmetry of the expansion coefficients, c^'"(j), under j 1 — J and the 
relation ()A.6|) between the Legendre functions one can rewrite Eq. ()4.38|) in the form 



= A-(x, + E 4(j>r (i) + E 1 dj cot .J 4-{j)vf' {g^ 

j — 2 Q, — di "^/^ 



(4.47) 



where 



v'fXx) = e{x - 1) {x" - l)"^/2p-_7(a;) . (4.48) 

Thus, one sees that the contribution of the integrals in Eq. ()4.38|1 at L = 1 in the region 
|x| < 1^1 is given by the term A^(x), which is determined completely by the values of the 
function at the points x = Obviously, the ffist two terms in ()4.47|1 have support in 
the ERBL region, while the integrals live in the DGLAP region. The term A^(x) can be 
restored from the knowledge of the function in the DGLAP region. 

At the scale /i' the function <f^{x) is given by the expression ()4.H8j) . To rewrite it in the 
form ()4.47|) one has to find the expansion coefficients c^(j) and c^'^(j) at the scale /i'. To 
this end one should insert the expansion (j4.38|) into (j4.39p and evaluate the corresponding 
integrals. This results in the following expression for the expansion coefficients at scale /j,' 



^.^U^^^l (4.49a) 

^ '^^Jcjj (j-j)(j+J-l) 



c 



c^'^if) + (-l)^-^c'^'-(/) 



^ Jci (j -j)(j +J-l) 

where the integration contour in both cases follows the line parallel to the imaginary axis 
such that 1 < Re j' < j (contour C/), and max{l, Rej, 1 — Re j} < Rej' (contour Cu). We 
notice that the integration contours can not be closed in the right half-plane because the 
integrals over the large semicircle do not vanish. Next, let us check that at /i' = (L = 1) 
the r.h.s reproduces the l.h.s. To this end let us shift the integration contour to the line 
Re j' = 1/2. Then the r.h.s of Eq. (j4.49a|) will be given by the sum of the residue at the point 
j' = j and the integral over the line Rej' = 1/2. The integral vanishes due to antisymmetry 
of the integrand under j 1 — j while the residue gives the coefficient c^'^(j). 

One sees that the coefficients c^(j), c^'^(j) mix under evolution [s^l- Therefore, their 
scale dependence is not independent. Nevertheless, it is straightforward to check that the 
transformation ^ = f{fi',c^) generated by Eqs. ()4.49|) possesses the necessary group 
property, /(/x", c'^) = /(/i", /(/.', c'^)). 

The coefficients (j) get contributions from c^'^(j). This effect can be interpreted as 
migration of partons from the DGLAP region to the ERBL region 0, E|- 

Further, let us note also that it follows from the Eq. (j4.49j) that for integer j > 2 

cf(j)+cf'+(j) + (-ir'cf'-(j)] = [4{j) + cl'\j) + i-iy-'c^^'-ij)]L-^^^K (4.50) 

It is not surprising that the combination in the brackets is nothing but the Gegenbauer 
moment of the GPD 

C^O = -V, I C^f. {£) ^lix) . (4.51) 
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FIG. 3: Example for the evolution of the function (f^=oA{x) = (1 — x^)'^''^ (curve (1)) and its 
derivative (curve (2)) to the scale L = 1.5, curves (3) and (4), respectively. 



For the physical GPD, H{x,C,), moments C.^C^jiO have to be polynomials in ^ of degree A^. 
This is the so-called polynomiality conditions It means that the knowledge of the GPD 
in the DGLAP region allows to restore the GPD in the central region, up to terms which 
satisfy the polynomiality conditions. 

Next, all singularities of the integrand in the integral in Eq. (|4.49a|) lie to the left of the 
integration line. This allows to shift it arbitrarily far to the right. Since the anomalous 
dimension 7(j) behaves as logj for large j we conclude that the coefficients c^'^(j) vanish 
faster than any power of 1/L for L oo. As far as the function v?g(x) in the DGLAP region 
|x| > 1^1 is expressed solely in terms of the coefficients c^'^(j) this implies also that the 
function in this region vanishes faster any power of 1/L for L —>■ oo. Thus the asymptotic 
expansion for the function (x) for ^' ^ oo (large L) has the form 



O . (4.52) 



Since the anomalous dimension 7(7) vanishes for j = 2 the function (x) tends to its 

asymptotic form, p^2^ ~ (1 — x^/^^), (see Refs. [H, 13), for 00. We stress that 

the expansion (j4.52|) is an asymptotic expansion and that the series in (I4.52p diverges for 
any x. 

Let us analyze now the behavior of the function ip'^ (x) at the point x = As follows 

from the properties of the functions qj^\x) the GPD (p'^ (x) is continuous at the points 
X = ±^ and has the following values (see Eqs. ()4.45|) and ()4.49a|) ) 
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where the k > 1. Taking into account the Eqs. (jA.ljl and ()A.7j) one finds that the first 
derivative of the function cp'^ (x) has at least a logarithmic singularity at x = even if the 
function ip^{x) is a smooth function at the input scale /i (see Fig. E)) 

^¥'f(x)^^^^a±|^rMog|x/|^|Tl| . (4.54) 

Here 

1 pK+ioo 

a± = -— dj'c^^^{j')L-<^'\ (4.55) 

If the integral in ()4.55|) diverges it means that the derivative of (p'^ (x) has a more singular 
behaviour at x = It is clear that the convergence properties of the integral in (j4.55|) is 
determined by the behavior of the coefficients c'^'^(j) for j — >■ ±2oo. The latter depends on 
the smoothness of the function ^p^ix) at the scale n, i.e. on the behavior of the function at 
the points x = and C, = ±1- 

The amplitude of a scattering process (DVCS, light meson production, etc.) is given by 
the convolution of a GPD with a hard scattering amplitude. The typical integral occurring 
in such a convolution at lowest order in as is 

/oo in^' (x) 

dx-^^^^\ , (4.56) 

where due to the support properties of the GPD the integration is restricted to the interval 
[— 1, 1]. We assume here that ^ > 0. Substituting (p^{x) by the expansion ()4.38|) and changing 
the order of summation (integration) over j and x one derives with the help of the Eqs. ()A.4|) 

j^JU-^) mJcSimrj j(j - 1) 

where the integration follows the line Rej = k, 1 < k < 2. 

Let us compare the approach presented here with the one developed in ji^ . One can see 
that the only difference between Eq. (81) in Ref. ji^ and our representation of the GPD, i.e. 
Eq. fl4.H8|) at L = 1, (or analoguosly Eq. (4) in Ref. [^) is, that the sum over discrete j in 
Eq. fl4.38|) is rewritten as an integral over j. To this end one has to construct the analytical 
continuation of the coefficients c^{j) with suitable analytic properties. This has been done 
in [i^. However, Eq. (69) in Ref. [i^, which provides such an analytical continuation, 
involves the function ip^{x) (or w{x, ,)Tn the notations of Ret. 0)ou«de the region where 
it is originally defined. Thus, such reformulation is possible only under certain additional 
assumptions on the analytic structure of GPDs. 



C. Small ^ asymptotics 

In this subsection we want to discuss the structure of the solution ()4.38|) in the limit 
^ — s> 0. (It is clear that when ^ ^ 1 one recovers the conventional ERBL evolution ^ 2o|.) 
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Obviously, the form of the function in the small ^ region can not differ strongly from the 
input function at L ~ 1. However, one may hope that the form of the evolved function for 
large scales, L ^ 1, depends weakly on the initial profile. 

The sum in Eq. ()4.38p contributes only to the central region. To estimate this sum one 
has to have some information on the expansion coefficients cg(j). If the coefficients c^{j) 
are not singular as ^ —^0, then at large L one can omit all the terms in sum except the 
first one. This type of behavior holds for some models of GPDs js^, 3^ but the opposite 



situation is not excluded. Thus, we restrict ourself to the analysis of the small ^ behavior 
of the integrals in Eq. ()4.38p . 

First of all we shall show, that in the limit — * and x fixed, Eq. (j4.38p takes the 
standard DGLAP form. Till the end of this section we shall assume that ^ > 0. Let us 
consider the integral in (|4.38|) corresponding to the term "a = +". For ^ — > and x fixed, 
the argument of the function gj^^ goes to infinity and so we can replace it by its asymptotic 
value 

The contour of the integration over j follows the line Rej = k, k > 1. The function gj^^ 
vanishes as for ^ and the parameter k can be chosen sufficiently large. In its turn 
the coefficient c^(j) becomes singular at ^ — 0. Indeed, using Eq. ()A.10|) one can derive 
the following representation 

where the integration goes from —ioo to +ioo staying to the right of all singularities of the 
integrand and 

!f^{s)= I dxx'-^^^{x). (4.60) 
Jo 

If Res is sufficiently large then (/^^(s) — * v^^=o(s) for ,^ — > even if the function v'g(x) 
is singular at a; = ^ = 0. The leading ^ asymptotics of c^(j) is determined by the 
rightmost singularity of the integrand in ()4.59p which is situated at s = j — 1. Thus one gets 

< w - - (I) ^i^^^o- - 1) (1 + ■ (^-61) 

Therefore, for ^ = the solution (jO^ {x > 0) takes the DGLAP form 

VIU^) = ^J/3 n=oU) L--'^'^'^ . (4.62) 

The integrals in ()4.38|) depend on the three parameters x, ^ and L. Eq. ()4.62|) was obtained 
in the limit ^ — 0, x and L fixed. Let us now keep the ratio x/C, = r > 1 fixed (i.e. we shall 
consider the DGLAP region) and evaluate the integral 

J(e,L,r) = - f djcl'+ij)L-'^^\f\r) (4.63) 
7r« Jc 
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for large 1/^ and L. The answer depends on the relative size of L and ^. Lets start with the 
situation when logL ^ log 1/^. In this case taking into account Eq. ()4.6H) one can evaluate 
the integral by the saddle point method. The position of the saddle point is determined 
from the equation 

1 - y (j;)cr = , a = log L/ log (2/0 . (4.64) 

For large a, (which of course corresponds to an unphysically high scale /x') = ACp/boa 
and the integral ()4.63p scales with L as 

J^L-^^^/^o(^°s>-i)g];^(r). 

Note that we restricted our consideration to the DGLAP region, x/^ > 1, only because 
the integrand in this region has no singularities in the half-plane Rej > 1, so that that 
we can freely deform the integration contour. In the ERBL region the integrand has poles, 
(see Eq. (j4.42p ). so when deforming contour one has to take into account the corresponding 
residues. These residues together with the terms in the sum in ()4.38|) result in the asymptotic 
expansion (|4.52j) . Consequently, the number of terms in the asymptotic expansion (|4.52p 
is controlled by the position of the saddle point ~ log L/ log (2/0- This means that 
the "time" r = logL which is necessary for the function to reach the asymptotic regime is 
proportional to log 2/^. 

Now let us decrease L or ^. Then a and decrease as well. When a <^ 1 (i.e. when the 
function is far from its asymptotic form) one finds 



+ (4.65) 
and recovers the well known "double scaling" behaviour 

log^/^(2/0 ' ' 

Let us now consider the situation when the function (p^=q{x) is enhanced at small x. 
Recent fits of parton densities (see Ref. Js^) shows that in the x region, 10~^ < a; < 10~^, 
they display a power-like behavior ip^=Q{x) ~ In this case the Mellin transform of the 
function (p^=o{x) has a pole at s = a. (In the case a > the first pole of (p^=o{j — 1) lays to 
the right of the first pole of the anomalous dimension 7(j).) At the same time the function 
(p^{x) is a regular function of x for ^ 7^ 0, so that the first pole of its Mellin transform lies 
at s = 0. It can be seen that near the point s = a the Mellin transform f^{s) takes the 
following form (p^{s) ~ (1 — ^''~"/(^, s))/(s — a), where /(^,a) = 1. One sees that the 
numerator contains two terms with different ^ dependence. Thus the integrand in ()4.63|1 
near the point j = 1 + a can be split up into two terms, which are proportional to 
and respectively. Each of these terms has poles, which, of course, cancel in the sum. 
Calculating the asymptotics of the first term (O""*) finds that in the case < 1 + a it 
is necessary to take into account the contribution due the pole at j = 1 + a which is 

I^°"=~r"gSa(^) (4.67) 
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For < a < 1 the anomalous dimension 7(1 + 0) is negative so that the evolution results 
in an enhancement of this contributions. The second term has the same ^ dependence as 
the pole contribution, However, this term vanishes faster with L than the first one. 

So one can expect that the pole contribution ()4.67|) will dominate in the limit L —>■ 00, 
logL/ log 1/^ fixed. Our result ()4.fi7j) agrees with the result of Shuvaev et al. (One can 
easily check that the integral (22) in Ref. [i^ gives rise to the function q^j^x^ (x/^).) 

V. SINGLET CASE 

In the singlet case one is dealing with the evolution of the quark, J-''^'~^{x,C,), and gluon, 
J^^{x,^), GPDs which are related to matrix elements of the C— even quark and gluon oper- 
ators, (see Eqs. ()2.5|) . ()2.25|) ). We shall proceed along the same lines as in the case of the 
isovector quark GPD. Similarly, we start with the formulation of the problem in coordinate 
space. The only difference with respect to the previous case is that we have now two func- 
tions f^izi, Z2), f^{zi, Z2) (see Eqs. ()2.24p ) instead of one {(p{zi, Z2)). For the sake of brevity 
we shall denote by /(-Zi, Z2) the two component function /(^i, Z2) = (/"^(-Zi, Z2), f^{zi, Z2)) = 
(f^^\zi,Z2),f^'^\zi,Z2))- Next, let us introduce a function ip^ip = {ip'^^ip^) = (^p^-^\ ijj^'^^) , as 
follows 

^P^'\z,,Z2) = (^i-Z2)V^'^(^i,^2), (5.1a) 

where k = 1,2. The functions ■ip'^^s) transform under SL{2, R) transformations according to 
Eq. (j4.3|) with s = 1/2 while the functions f^^^^ obey the same transformation with s = 1 
and s = 3/2, respectively. We define the scalar product on the space of functions / as follows 

(/1I/2) = {Hi^2) = + , (5.2) 

where the scalar product ('^/'f IV'?) the standard scalar product on the space L^(i? x R) , 
(see Eq. (gH)). 

Since the functions -ipi^^^ obey the SL{2, R) transformations ()4.3p with conformal spin 
s = 1/2 they can be expanded in eigenfunctions of the Casimir operator according to the 
Eq. ()4.21|) . Again, we shall omit the Fourier integral over ^ and write down all expansions 
for the function f^{z) (ip^^^z)) 

f{zu Z2) = I ^e-«(^^+^^) f^{z, - Z2) . (5.3) 

Obviously, an expansion for the function ip^lx) has the form ()4.2(i|l . Then we shift the 
integration contour from the line Rej = 1/2 to Rej = 5/2. Doing so one has to take into 
account the residues at the points j = 1, 2 which cancel identically with the first two terms in 
the sum over j. To see that the sum of the residues at j = 2 for the function ip^{z) vanishes 
one should remember that this function is symmetric with respect to 2; — > —z. Thus the 
expansion for the function f^{z) can be written in the following form 

f^'\z) = J2^ij)4\j)^'/"\^) - I / £-M3)af^\3)<^f'\z^), (5.4) 
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where 2 < k < 3, k = 1,2, z± = ±z + iO and the coefficients a^^\j),a^^^''^{j) are given by 
Eq. ()4.23|) after the substitution tjj — > ip'^'^^ The functions '^f^^\z±) are defined as follows 



(and similar for '^^''^^\z)) where ^E'|(-2;) is defined in Eq. ()4.15|1 . Now, using the same 
arguments as in Sect. IIVI one can show that the Hamiltonians describing the renormal- 
ization of the singlet operators Eqs. ()2.17|) - ()2.21|) act on the functions (^f\zi,Z2) = 
e-i?(^i+22)^|.(^0(^^^^) follows 

W'<^f\z,,Z2) = E^\J)<^f{z^,Z2). (5.5) 

Note that in the formula above no summation over repeated indices i, k is implied. The 
matrix of the "anomalous dimensions" E"^^ has the well known form 

E^'^ij) ^ E'\j) =Cp {2 [^(j) - ^(2)] - -^-^ + , (5.6a) 

^^^^^^-^^^^^^-^ (..4-l)(.-2) - ^^•^^) 
E^\j) ^ E'\j) =c /.';.^\'^ , (5.6c) 



E^\3) = E'\3) =2iV, 



V'(j)-^(3)- ^ ^ 



\Nc^\ns. (5.6d) 



j(j + l) (j-l)(j-2) 

The solution of the evolution equation (j2.16j) for the function /^(-z) can be written in the 
form 

00 



i=3 



- I [ -^^(J) {L-^'^H'^U)), ^f\z±) ■ (5.7) 
2 sm vrj * ' k j 

Here, 7(j) is the matrix of the anomalous dimensions, 7(7) = 2E{j)/bQ, a^(j) is a two 



dimensional vector, a't{j) = {a'^'^^\j) , ac'^'^\j)) , and the integration contour is the same as 



in Eq. 

By taking the Fourier transform Eq. (j5.7p can be brought into the form 

{) = E ft (lii ) i-"'' ^n^) + E X I (i) <^ro) ■ (5-8) 

The quark and gluon GPDs J-'^^-^i^^O defined in Eq. ()2.26p . The matrix functions Pj{x) 
and qj{x) are given by the following expressions 

Pj(^) = I n ^ (2)/ A ) (5-9) 



epr(^) 

^\x) 
egf(x) 
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and the functions p^™'\qj^^ are defined by Eqs. (j4.4m) . (j4.4H) . The expansion coefficients, 
c= (c'-^^c'-^^) have the following form 



dxPi-i { — 



X 



dx'^ " 
d^ 



(5.11a) 



(5.11b) 



Note, that the gluon GPD J^2{x,^) = J^g{x,^) is a symmetric function of x, while the 
quark GPD J^i{x,^) = J-'q{x,^) — J^g{—x,^) is an antisymmetric one. Therefore, all integer 
moments with even j (Eq. (j5.11a|) ) vanish and c'^'~^{j) = c^'~{j). The matrix exponent 
L^'^^^^ can be written in the form 



where 



(5.12) 



(5.13) 



A(j) 



(711 (j) - 722(j))^ 



+ 7l2(j)72l(j) • 



(5.14) 



The projectors A±{j), A\{j) = A±(j), A^{j)A^{j) = are given by the following expression 



I=F 



711 (j) -722 (i) 
2 



A(j) \ 721 (j) 



712 (j) 
722(i)-7ii(i) 



(5.15) 



Let us note that although each term on the r.h.s of ()5.12|1 has branching points, their sum 
is an analytic function. For large j, 7+(j) ~ 711 (j) = 7,(7) and 7_(j) ~ 722(7) = Igij) so 
we shall refer to 74.(j)(7_(j)) as the quark (gluon) anomalous dimension. Let us notice that 
due to mixing the quark (but not gluon) anomalous dimension, 7+(j), has a pole at j = 2. 

The formulae ()5.8|1 . 1)5.111) represent the solution of the evolution equation for the quark 
and gluon GPDs. 

Having set L = 1 in Eq. ()5.8|) one derives a representation analogous to that in the 
isovector case, Eq. ()4.47|1 . 



^^(x, = A^{x, + p^P, c-^(j) + Y.jij^ cot TTj v; 



ax 



c^^'^ij), (5.16) 



where the integration follows the line Rej = 1/2. The function Vj is defined as 





±P)'\x) 





^rl'\x) 



where 



Vp{x) = 9{x - 1) (x^ - l)^/2p-_fe , 
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x] 



(5.17) 



(5.18) 



The vector A^{x,C,) is given by the following expression 



X 



A^(^,0 =^(iei-kl)fe(ieU) + 



[X 



m 



dx ^ 



(leu) 



(5.19) 
(5.20) 



The expansion coefficients at the scale /x' are given by the same expressions as in the non- 
singlet case, Eqn. ()4.49|) . which should now be understood as vector equations. 

The anomalous dimension 7+(jQ vanishes at j = 3, so that for ^ oo the GPD takes 
the well known asymptotic form p, 3^ 



X 



■as 
2 



c 



1 



X 



Uf— \l- — \D 



2N 2 



1 - 



X 



D 



where 



D 
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(5.21) 
(5.22) 

4(n, + 4Cp)7.r^"^"^'^^ ■ -^v-s.. (5.23) 

Similarly to the case of the isovector GPD one can analyze the small ^ behavior of the 
quark and gluon GPDs. In the situation when the GPDs are the smooth functions at the 
input scale their form after evolution in the small ^ region is determined by the rightmost 
pole of the anomalous dimension 7+(j), which is located at j = 2. The corresponding 
asymptotics have the "double scaling" form which results from evaluation of the integral in 
Eq. ()5.8|) by the saddle point method. We again consider the integral corresponding to the 
term with a = + in Eq. (j5.8|l . The expansion coefficients (j5.1H) can be expressed as follows 



J dx[xJ^i{x,^) + J^2ix,^)] . 



m{j) 



m{j) 



r(j + 1/2) 

0Fr(j-i) 



--^l(j-l) 

U - 2) J-2(j - 2) 



(5.24) 
(5.25) 



where 



f{s)= [ dxx'-^f^'{x,0 ■ (5.26) 
Jo 

Inserting ()5.24|1 into Eq. ()5.8j) one finds that the saddle point is determined by the equation 

logL 



1 -7-(j*)f^ 
For small a one finds j = 2 + e with 



iog2/e 



(5.27) 



(5.28) 
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Then one finds for tlie integral X (wliicfi corresponds to the term a = + in Eq. ()5.8p ) in the 
"double scaling" approximation 

^) " (^) 2 .^^^^^ ^^^2 + e) m(2 + e). (5.29) 

Here r = x/C, and we assume that ^ > 0. Due to symmetry properties taking into account 
the term a = — results in the replacement 

in the above expression. 



VI. SUMMARY 

We developed a method for solving the leading order evolution equations for Generalized 
Parton Distributions. The form of the solution, Eqs. ()4.38|) and ()5.8|) . is fixed completely by 
the symmetry properties of the evolution kernels. We have shown that the GPD should be 
treated differently in the ERBL and DGLAP regions, which reflects both the mathematical 
structure of the equation and its physical content. 

The Eqs. (I4.38j) and (j5.8j) can be used for both analytical and numerical studies of the 
GPD's evolution. It is important to stress here that all quantities entering these equations 
are unambiguously defined. We demonstrated that in the limits ^ — > 1 and ^ — > the so- 
lutions ()4.38|) and ()5.8|) of the GPD evolution equations take the form of the solutions of 
the ERBL and DGLAP evolution equations, respectively. Furthermore, using the represen- 
tations ()4.38p and ()5.8p one can easily obtain the asymptotic form of the GPDs in different 
regimes (large L, small ^, and so on). 

It would be quite interesting to apply the developed approach to the analysis of the 
analytic structure of the anomalous dimensions of twist-3 operators. It is well known ||4j| 
that the knowledge of the energies of the reggeon bound states allows to extract anomalous 
dimensions of the operators at the singular point of in the j— plane. This correspondence has 
been checked for the twist two-operators. Recently, methods for the calculation of multi- 
reggeon bound states were developed ji^ EE and the predictions for the anomalous 



dimensions of higher twist operators at the singular points were obtained |47|, |48[ . However, 
so far it is not known how to solve the problem of the analytical continuation of the anoma- 
lous dimensions of higher twist operators. We hope that our approach can be generalized to 
higher twist operators, at least for the class of the twist-3 operators for which the evolution 
equation are known to be integrable 27, 2^ IE 5^ - 
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APPENDIX A: AUXILIARY FORMULAE 



In this Appendix we collected formulae which were useful in our analysis. We use the 
definition of the Bessel and Legendre functions of Ref. It follows from the properties 
of the Legendre functions that the functions p^p\x), qf^\x) defined in ()4.4()j) . ()4.41|1 satisfy 
the following relations 

^pT\^) = pT'"\^). (A. la) 

= qf'~'\x). (A. lb) 

The function qf^\x) can be represented in the form 

^-\x) y^'QT\^+^^)-^-'^'QT\^-'^) ^ (A.2) 

^ 2i sin irj ' 

where 

Qf\z) ={z' - ir/' Q-!^,{z) . (A.3) 
We give here some integrals involving the functions pj™""* (x) and gj'"'* (x) 



pf\x) = {-l)^2'^T{m)^. — ^, (A.4a) 



il — X ■' r(j + m 

dx _ 2™r(m)r(j -m) 



q'^Ux) = (A.4b) 

1 x + ^ ' 2sin7rj r(j + m) ^ ' 

and 

f , J r\^) = (-l)"^2-r(m)-^I|^^ . (A.4c) 
_i X — l±%\j 1 sm TTj 1 (j + m) 

The following integrals involving the Legendre functions were useful 

, ^ / X ^ / ^ 2 (-l)"^sin7rz/ ^ , 

dxP^ix) Pm{x) = —. — (. ^ — , A.5a 

OO 1 

dxP,{x)Q,{x) = (A.5b) 

I (z/ - A)(z/ + A + 1) 

where m is integer. Next, we remind the following relation 

P,{z) = -tanvrz/ [Q,{z) - Q-,-i{z)] (A.6) 

71 

between the Legendre functions of the first and second kind. Next, for z ^ 1, the Legendre 
function Qu{z) has the following behavior 

Q^{z) = -hog{z-l) + 0{l). (A.7) 
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It follows from the Eqs. ()A.5b|) and ()A.6j) that 

dxP_i/2+ipix)P^i/2+ip'{x) = cothTrp [S{p - p') + 5{p + p')] . (A. 8) 

The Fourier transform of the function "^^jiz) (see Eq. fl4.15|) ) is 

where Rej > m. We also remind that 



(A.9) 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7 

[8] 
[9] 

[lo; 
[11 

[12 

[13 

[14 
[15 

[16 
[17 

[is: 

[19 

[20 
[21 



dxx-^p^{x) = r-' ^ 'I ) ' ^ , (A.IO) 

TT^/'^L (S) 







dxx ^P^{x)-2'' Y ■ (A-11) 



A. Manashov, M. Kirch and A. Schafer, Phys. Rev. Lett. 95, 012002 (2005) 

D. Midler, D. Robaschik, B. Geyer, F. M. Dittes and J. Horejsi, Fortsch. Phys. 42, 101 (1994) 

X. D. Ji, Phys. Rev. Lett. 78, 610 (1997) 

A. V. Radyushkin, Phys. Lett. B 380, 417 (1996) 

M. Diehl, Phys. Rept. 388, 41 (2003) 

K. Goeke, M. V. Polyakov and M. Vanderhaeghen, Prog. Part. Nucl. Phys. 47, 401 (2001) 
A. V. Belitsky and A. V. Radyushkin, "Unraveling hadron structure with generalized parton 
distributions," |arXiv:hep-ph /0504030| 



A. Airapetian et al. [HERMES Collaboration], Phys. Rev. Lett. 87, 182001 (2001) 
S. Stepanyan et al. [CLAS Collaboration], Phys. Rev. Lett. 87, 182002 (2001) 
A. Ahapetian et al. [HERMES Collaboration], Phys. Rev. Lett. 94, 012002 (2005) 
P. Hagler, J. Negele, D. B. Renner, W. Schroers, T. Lippert and K. Schilling [LHPC collabo- 
ration], Phys. Rev. D 68, 034505 (2003) 

M. Gockeler, R. Horsley, D. Pleiter, P. E. L. Rakow, A. Schafer, G. Schierholz and W. Schroers 
[QCDSF Collaboration], Phys. Rev. Lett. 92, 042002 (2004) 

P. Hagler, J. W. Negele, D. B. Renner, W. Schroers, T. Lippert and K. Schilling [LHPC 

Collaboration], Phys. Rev. Lett. 93, 112001 (2004) 

M. Gockeler et al. [QCDSF Collaboration], arXiv:hep-lat/0501029l 

A. V. Radyushkin, Phys. Lett. B 385 (1996) 333; I. I. Balitsky, A. V. Radyushkin, Phys. Lett. 
B 413 (1997) 114. 

V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972) [Yad. Fiz. 15, 781 (1972)]. 

Y. L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977) [Zh. Eksp. Teor. Fiz. 73, 1216 (1977)]. 

G. Aharelh and G. Parisi, Nucl. Phys. B 126, 298 (1977). 

A. V. Efremov and A. V. Radyushkin, Phys. Lett. B 94, 245 (1980). 

G. P. Lepage and S. J. Brodsky, Phys. Rev. D 22, 2157 (1980). 

N. Kivel and L. Mankiewicz, Nucl. Phys. B 557, 271 (1999) 



33 



[22] A. Shuvaev, Phys. Rev. D 60, 116005 (1999) 
[23] J. D. Noritzsch, Phys. Rev. D 62, 054015 (2000) 

[24] A. V. Belitsky, B. Geyer, D. Miiller and A. Schafer, Phys. Lett. B 421, 312 (1998) 

[25] D. Miiller and A. Schafer, arXiv :hep-ph/0509204 

[26] I. I. Bahtsky and V. M. Braun, Nucl. Phys. B 311, 541 (1989). 

[27] V. M. Braun, S. E. Derkachov and A. N. Manashov, Phys. Rev. Lett. 81, 2020 (1998) 
[28] V. M. Braun, S. E. Derkachov, G. P. Korchemsky and A. N. Manashov, Nucl. Phys. B 553, 
355 (1999) 

[29] S. J. Brodsky, Y. Frishman, G. P. Lepage and C. T. Sachrajda, Phys. Lett. B 91, 239 (1980). 

[30] L. Mankiewicz, G. Piller and T. Weigl, Eur. Phys. J. C 5, 119 (1998) 

[31] V. M. Braun, G. P. Korchemsky and D. Miiller, Prog. Part. Nucl. Phys. 51, 311 (2003) 

[32] I.M. Gelfand, M.I. Graev and N.Ya. Vilenkin, Generalized functions. Vol.5 Integral geometry 

and representation theory, New York, NY Academic Press 1966. 
[33] Vilenkin, N. Ja. Special functions and the theory of group representations. Translated from 

the Russian by V. N. Singh. Translations of Mathematical Monographs, Vol. 22 American 

Mathematical Society, Providence, R. I. 1968 
[34] F. Yuan, Phys. Rev. D 69, 051501 (2004) 
[35] A. V. Radyushkin, Phys. Lett. B 449, 81 (1999) 
[36] L. Pukanszky, Trans. Amer. Math. Soc. 100 (1961) 116-152 

[37] I. S. Gradshtejn and I. M. Ryzhik, Tables of Integrals, Sums, Series and Products, Academic 
Press, NY, 1980 

[38] I. V. Musatov and A. V. Radyushkin, Phys. Rev. D 61, 074027 (2000) 

[39] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 0207, 
012 (2002) 

[40] A. De Rujula, S. L. Glashow, H. D. Politzer, S. B. Treiman, F. Wilczek and A. Zee, Phys. 

Rev. D 10, 1649 (1974). 
[41] R. D. Ball and S. Forte, Phys. Lett. B 335, 77 (1994) Phys. Lett. B 336, 77 (1994) 
[42] A. G. Shuvaev, K. J. Golec-Biernat, A. D. Martin and M. G. Ryskin, Phys. Rev. D 60, 014015 

(1999) 

[43] T. Jaroszewicz, Phys. Lett. B 116, 291 (1982). 

[44] S. E. Derkachov, G. P. Korchemsky and A. N. Manashov, Nucl. Phys. B 617, 375 (2001) 
[45] S. E. Derkachov, G. P. Korchemsky, J. Kotanski and A. N. Manashov, Nucl. Phys. B 645, 
237 (2002) 

[46] H. J. De Vega and L. N. Lipatov, Phys. Rev. D 64, 114019 (2001) 

[47] G. P. Korchemsky, J. Kotanski and A. N. Manashov, Phys. Lett. B 583, 121 (2004) 

[48] H. J. de Vega and L. N. Lipatov, Phys. Rev. D 66, 074013 (2002) 

[49] S. E. Derkachov, G. P. Korchemsky and A. N. Manashov, Nucl. Phys. B 566, 203 (2000) 
[50] A.V. Belitsky, Phys. Lett. B 453, 59 (1999); Nucl. Phys. B 558 (1999) 259; Nucl. Phys. B 
574 (2000) 407. 

[51] The coefficients c'i'^{j) are the matrix elements of the so-called string operators [2^ . 



34 



